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The  general  theme  of  our  research  has  been  to  develop  stochastic  net¬ 
work  processes  for  modeling  the  movement  of  discrete  units  in  networks. 
Primary  examples  are  the  movement  of  parts  and  supplies  in  manufactur¬ 
ing  plants  and  in  distribution  systems  and  the  movement  of  data  packets 
and  telephone  calls  in  computer  and  telecommunications  networks.  The 
feature  of  our  research  is  the  emphasis  on  the  next  genera¬ 
tion  of  intelligent  networks  that  will  be  the  backbone  of  our  mamrfacturing 
and  computer  systems.  In  these  networks,  the  processing  of  units  at  the 
nodes  and  the  routing  of  units  typically  depend  dynamically  on  the  actu^ 
network  congestion,  and  units  move  concurrently  (e.g.  batA  processing). 
Most  of  the  present  theory  of  stochastic  network  processes  is  for  munte^- 
gent  networks  in  which  the  nodes  operate  independently,  the  routes  of  units 
are  independent,  and  the  imits  move  one-at-a-time.  Our  goal  is  to  provide 
an  understanding  of  these  more  complex  intelligent  networks  by  describing 
their  stochastic  behavior. 

Another  focus  of  our  network  research  is  on  parallel  simijlation,  which 
is  one  of  the  most  promising  areas  for  the  use  of  parallel  or  distribut^  pro¬ 
cessing.  One  task  is  to  develop  stochastic  network  models  for  assessing  the 
feasibility  and  quality  of  various  protocols  in  parallel  simulations.  Empha¬ 
sis  is  on  an  optimistic  or  aggressive  Time  Warp  protocol  that  periodically 
generates  superfluous  data  that  has  to  be  rolled  back  or  discarded.  Another 
task  is  to  develop  algorithms  that  can  be  incorporated  as  subroutines  in 
certain  types  of  parallel  processing  simulations,  such  as  queueing  ne^orks, 
in  which  the  system  evolution  can  be  represented  by  recursive  equations. 

The  following  is  a  summary  of  the  papers  we  have  written  for  this  project 
during  the  last  two  years. 


1 


Dist 


Availability  Codes 


Avail  and/or 
Special 


f 


1  Relating  the  Waiting  Times  and  Queue  lengths  in  Heavy  IVafiic 
Systems  by  W.  Szczotka,  W.  Topolski  and  the  PL  Stock.  Processes  Appl. 
52,  119-134,  1994. 

Diffusion  models  have  been  instrumental  in  analyzing  queueing  systems  in 
heavy  traflac.  We  show  that  there  is  more  than  the  one  conventional  reflected 
Brownian  motion  model  for  such  analyses  —  there  are  three  natural  models 
(including  the  conventional  one)  depending  on  the  subtle  nature  of  the  heavy 
traffic.  These  are  estabUshed  via  weak  convergence  limit  theorems.  Other 
major  results  show  that,  in  heavy  traffic,  an  arriving  umt  s  waiting  tune 
before  its  service  begins  is  approximately  equal  to  the  number  of  units  in  the 
system  the  average  service  time  (this  is  not  true  in  light  or  moderate 
traffic).  This  property  is  useful  for  obtaining  confidence  intervals  for  queue 
lengths  based  on  waiting  tune  mformation  or  vice  versa. 

2  Campbell’s  Formula  and  Applications  to  Queueing  by  Schmidt, V. 
and  the  PI.  To  appear  in  Frontiers  in  Queueing  1995. 

Campbell’s  formula  for  Palm  probabilities  is  a  basic  tool  for  deriving  prop¬ 
erties  of  stationary  queueing  systems  and  stationary  processes  in  general. 
This  study  shows  that  Campbell’s  formula  is  essentially  Fubini’s  theorem 
for  random  measures  —  this  clarifies  many  of  its  previous  ad  hoc  applicar 
tions  and  ccmjectures  about  its  scope.  We  also  give  new  insights  into  the 
versatility  of  Campbell’s  formula  by  establishing  several  equivalent  versions 
of  it.  A  few  of  these  are  known  (e.g.  the  exchange  formula  and  H  =  AG). 
Also  included  are  applications  involving  integrals  with  respect  to  random 
product-measures,  waiting  times  in  systems,  rate  conservation  laws,  sojourn 
times  of  processes,  travel  times  in  networks,  ladder  heights  in  risk  processes 
and  virtual  delays  in  queueing  systems. 

3  Performance  Limitations  of  Parallel  Simulations  by  Chen,  L.  and 
the  PI.  Submitted  for  publication  1994. 

This  study  shows  how  the  performance  of  a  parallel  simulation  may  be  af¬ 
fected  by  the  structure  of  the  system  being  simulated.  We  consider  a  wide 
class  of  “linearly  synchronous”  simulations  consisting  of  asynchronous  and 
synchronous  parallel  simulations  (or  other  distributed-processing  systems), 
with  conservative  or  optimistic  protocols,  in  which  the  differences  in  the 
virtual  times  of  the  logical  processes  being  simulated  in  real  time  t  are  of 
the  order  o{t)  as  t  tends  to  infinity.  Using  a  random  time  transformation 


2 


idea,  we  show  how  a  simulation’s  processing  rate  in  real  time  is  related  to 
the  throughput  rates  in  virtual  time  of  the  system  being  simulated.  This 
relation  is  the  basis  for  establishing  upper  bounds  on  simulation  processing 
rates.  The  boimds  for  the  rates  are  ti^t  and  are  close  to  the  actual  rates 
as  numerical  experiments  indicate.  We  use  the  boimds  to  determine  the 
maximum  number  of  processors  that  a  simulation  can  effectively  use.  The 
bounds  also  give  insight  into  efficient  assignment  of  processors  to  the  logical 
processes  in  a  simulation. 

4  Parallel  Simulation  by  Multi-Instruction  Longest  Path  Algo¬ 
rithms  by  Chen,  L.  and  the  PI.  Submitted  for  publication,  1994. 

This  paper  presents  several  basic  algorithms  for  the  parallel  simulation  of 
G/G/1  queueing  systems  and  certain  networks  of  such  systems.  The  cover¬ 
age  includes  systems  subject  to  manufacturing  or  commumcation  blocking, 
or  to  loss  of  customers  due  to  capacity  constraints.  The  key  idea  is  that 
the  customer  departure  times  are  represented  by  longest-path  distances  in 
directed  graphs  instead  of  by  the  usual  recursive  equations.  This  repre¬ 
sentation  leads  to  algorithms  with  a  high  degree  of  parallelism  that  can 
be  implemented  on  parallel  computers  with  single  or  multiple  instruction 
streams. 

5  Little  Laws  for  Waiting  Times  and  UtUity  Processes  by  the  PI. 
Queueing  Systems  17,  137-181,  1994. 

A  fundamental  law  that  holds  for  many  service  systems  is  that  the  average 
time  W  that  a  unit  waits  in  the  system  is  directly  proportional  to  the  average 
queue  length  L  of  the  system,  i.e.,  L  =  XW ,  where  A  is  the  average  arrival 
rate  of  units  to  the  system.  This  law  is  known  to  be  imiversal  in  the  sense 
that  it  is  true  when  all  three  terms  exist.  Establishing  the  existence  of  these 
limits  for  a  new  system,  however,  may  pose  a  problem  and  this  is  where 
there  are  stiU  open  questions.  This  study  gives  a  new,  comprehensive  and 
simpler  aproach  for  proving  Little  laws  (i.e.  establishing  the  existence  of 
the  limits).  This  includes  laws  for  general  utility  processes  as  well  as  for 
waiting  times.  Using  this  approach,  we  obtain  fundamental  Little  laws  for 
Markovian,  regenerative  and  stationary  systems  that  apply  to  large  classes 
of  systems;  prior  to  this,  there  were  only  ad  hoc  resffits  in  this  regard  that 
were  unfriendly  to  non-specialists. 
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6  Parallel-Processing  Times;  Extreme  Values  of  Phase-type  and 
Mixed  Random  Variables  by  Sungyeol  Kang  and  the  PI.  Submitted  for 
publication,  1994. 

In  computer  or  telecommunications  systems,  a  typical  concern  is  the  time  it 
takes  for  a  group  of  data  packets  that  will  eventually  constitute  one  message 
to  be  processed  by  a  network  of  computers.  This  time  is  the  maximiun  or 
cxtrtmt  voiu6  of  the  travel  times  of  the  units  through  the  network.  A  simi¬ 
lar  concern  is  the  time  to  complete  a  job  in  a  manufacturing  (or  computer) 
system  consisting  of  many  tasks  performed  in  parallel  by  a  network  of  work¬ 
stations.  An  example  is  that  20  units  must  complete  a  PERT  network  in 
parallel  before  they  are  brought  together  as  one  system.  A  basic  problem 
is  to  determine  the  distribution  of  the  time  to  complete  a  large  number  of 
tasks  in  parallel.  That  is,  knowing  the  probability  structure  of  the  individ¬ 
ual  tasks,  what  is  the  asymptotic  distribution  of  the  maximum  of  the  task 
t.iTTiPs  as  the  number  of  tasks  tends  to  infinity? 

We  address  this  parallel-processing  (or  extreme-value)  problem  for  the 
following  settings:  (a)  A  task  consists  of  performing  a  set  of  randomly  se¬ 
lected  subtasks  in  series  and  the  subtask  durations  have  Erlang  distributions, 
(b)  The  task  times  are  independent,  identically  distributed  phase-type  ran¬ 
dom  variables  (e.g.  completion  times  of  Markovian  PERT  networks  or  task- 
graphs).  (c)  The  tasks  are  dependent  and  their  distributions  are  selected  by 
a  random  environment  process  such  as  a  Markovian  or  stationary  selection 
process. 

We  show  that  the  distributions  of  the  task  completion  times  in  these 
settings  are  the  classical  extreme-value  distributions,  even  though  the  task 
timps  are  dependent  and  may  involve  a  series  of  dependent  subtasks.  In¬ 
cluded  are  three  fundamental  results  that  apply  to  general  extreme-value 
contexts: 

(1)  All  Phase-type  distributions  are  in  the  domain  of  attraction  of  the  Gum- 
bel  distribution. 

(2)  Extreme  values  of  Markovian-  or  stationary-selected  variables  are  deter¬ 
mined  by  certain  states  of  the  selection  process  whose  associated  distribu¬ 
tions  have  a  “dominant  tail” . 

(3)  Extreme  values  of  sums  of  random  variables  are  determined  by  parts  of 
the  sum  that  dominate  the  other  parts. 
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7  Partition-Balanced  Markov  Processes  by  Alexopoulos,  C.,  El-Tannir, 
A.  and  the  PL  Submitted  for  publication,  1994. 

When  can  the  stationary  distribution  of  a  Markov  process  be  obtained  by 
pasting  together  several  stationary  distributions  that  represent  the  process 
restricted  to  certain  subspaces?  This  study  describes  a  class  of  “partition- 
balanced”  Markov  process  that  have  this  cut-and-paste  or  divide-and-conquer 
property.  The  importance  of  this  property  is  that  the  problem  of  obtaining 
a  stationary  distribution  on  a  large  space  (e.g.  for  networks)  reduces  to 
finding  several  stationary  distributions  on  smaller  subspaces,  possibly  even 
by  simulations. 

The  notion  of  partition-balance  is  a  “macro-reversibility”  property  re¬ 
sembling  the  detailed  balance  property  of  reversible  processes.  We  present 
several  characterizations  of  partition-balance  and  identify  subclasses  of  tree¬ 
like,  starlike  and  circular  partition-balanced  processes.  A  new  “circular 
birth-death  process”  is  also  presented  as  a  vehicle  for  the  analysis.  The 
results  are  illustrated  by  a  queueing  model  with  controlled  service  rate,  a 
multi-type  service  system  with  blocking  and  a  parallel-processing  model.  A 
few  comments  address  partition-balsmce  for  non-Markovian  processes. 
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